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ABSTRACT 

The formation and evolution of massive red galaxies form a crucial test of theories of galaxy formation based 
on hierarchical assembly. In this letter we use observations of the clustering of luminous red galaxies from the 
Bootes field and N-body simulations to argue that about 1 /3 of the most luminous satellite galaxies appear to 
undergo merging or disruption within massive halos between z — 0.9 and z — 0.5. 
Subject headings: galaxies: evolution - galaxies: halos 
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1. INTRODUCTION 

The assembly of the most massive galaxies is a key test 
of cold dark matter (CDM) models of galaxy formation, as 
the ongoing growth of massive galaxies via mergers is a 
generic feature of hierarchical CDM models. Observation- 
ally, the most massive galaxies have little ongoing star for- 
mation, an d the bulk of their stellar mass was formed at 
z > 2 (e.g. iBower. Lucev & Ellisl Il99l iTrager et all 120001: 
ICool et al.l 120061 and references therein). If there is appre- 
ciable growth of these galaxies at z < 1, this must be due to 
galaxy mergers, as predicted by the hierarchical CDM mod- 
els. 

Evidence for the ongoing assembly of massive galax- 
ies is inconclusive. While the stellar mas s within the red 
galaxy population has doubled since z = 1 dBell et alJl2004l 
IWillmeretaT]|2006t iBrown et al1l2006al) . this appears to be 
due to the truncation of star format ion in blue galaxies, and 
the role of merger s is poorly known. Ivan Dok kum ( 2005b and 
iBell et all (l2006allbh . using close galaxy pairs, conclude that 
L+ red galaxies grow rapidly via mergers since z = 1, while 
Masjedietal. (2006), using similar techniques, find that the 
merger rate of 4L* red galaxies is only ~ 1% Gyr -1 . Using the 
galaxy space density, IBrown et all (l2006al) find that the stel- 
lar masses of 4L* red galaxies grow by ~ 25% since z ~ 0.7, 
while others find no significant g r owth over similar redshift 
ranges (e.g.. iBundv et alj 120061: ICimatti. Daddi & Renzinil 
2006; iCaputi et alJl2006tlWake et alJl2006l) . 

There is an additional route to constraining the evolution of 
galaxies, which is to use their clustering properties. Building 
upon the theoretically understood evolution of the dark mat- 
ter halo population, we can obtain complementary constraints 
which bypass the model dependence of stellar evolution or 
merger times as a function of projected distance. We illustrate 
this approach in this Letter, presenting evidence from the evo- 
lution of their clustering that luminous red galaxies undergo 
merging or disruption between z ~ 0.9 and z ~ 0.5. 

2. THE OBSERVATIONAL SAMPLE 

We use galaxies in the 9 deg 2 Bootes field, which has been 
imaged in the o ptical and infrared b y the NOAO Deep Wide- 
Field (NPW FSj Jannuzi & Devil 99% and Spitzer IRAC Shal- 
low Surveys dEisenhardt et all 20041) . We use a luminous (> 



1.6L+) subset of the Bootes red galaxy sample which was se- 
lected from the imaging using empirical photomet ric redshifts 
and an evolving restframe U — V color criterion (IBrown et all 
2006a). This subset includes galaxy samples in three redshift 
slices: 0.4 < z < 0.6, 0.6 < z < 0.8 and 0.8 < z < 1 .0 with co- 
moving volumes of 2.4, 3.5 and 5.2 x 10 6 (/i _1 Mpc) 3 respec- 
tively, and at each redshift, galaxies are selected to be above 
a luminosity threshold such that the sample has a constant 
comoving number density (n = 10~ 3 /z 3 Mpc~ 3 ). Our results 
are based on the observed evolution of the angular clustering 
of these samples, containing a few thousand galaxies each. 
We transform from models of the spatial clustering to angular 
clustering using a redshift distribution model which accounts 
for the small measured uncertainties of the photometric red- 
shifts (<r z < 0.05). We de scribe the clustering measurements 
and modeling in detail in lBrown~e t al. ( 20063. 

3. MODELING GALAXY CLUSTERING 

Our galaxy samples have been selected to have constant 
h. We find that their clustering evolves very little, with 
£(6/!~ 1 Mpc,z) ~ 1. We begin with an analytic argument to 
show that these galaxies cannot be undergoing a common lu- 
minosity evolution history, such as pure passive evolution, 
without any mergers 1 . In such a scenario each galaxy pre- 
serves its identity and no galaxy leaves or enters the sam- 
ple. If we assume that galaxies and mass follow the same 
velocity field, and retain their identities, then the continuity 
equat ion in the linear regime requires that J ga i = i5 m dPeeblesI 
1980]). If we define S g . d \(z) = b(z)S m (z) and the growth func- 
tion/^) ee S m (z)/S m (0) then b(z) = 1+Er\z) [b(Q)-l] (H3 
1996). As shown in Figure[T] this prediction is in good agree- 
ment with our numerical simulations with passive evolution 
(see § 13.31 ). Assuming scale-independent, deterministic bias- 
ing therefore predicts evolution in £ which is not in agreement 
with the observations. In fact, we find that passive evolu- 
tion cannot fit the trend of the central values of the clustering 
strength for any plausible cosmology. 

To go further we need a way of connecting the galaxies we 
observe with the host dark matter halos whose evolution the- 

1 M.W., M.B. and A.D. thank Ravi Sheth for discussions at the Aspen 
Center for Physics which generated this argument. 
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FIG. 1 . — (Top) The correl ation length, rp, from a power-law fit for the 
sample as a function of z from Brown et al. (2006b). (Bottom) The evolution 
of the large-scale bias, assuming a dark matter power spectrum with <r 8 = 
0.8. The solid line is b(z) = 1 +£>"' [b(0) — 1] (see text), the open circles are 
measured from our "passive mocks" and the solid circles are from the mocks 
which best fit w(8) at each redshift. 

ory predicts. The halo model (see e.g. Cooray & Shethl 120021 
for a review) has provided us with such a physically informa- 
tive and flexible means of describing galaxy bias. The key 
insight is that an accurate prediction of galaxy clustering re- 
quires a knowledge of the occupation distribution of objects in 
halos (the HOD) and their spatial distribution. In combination 
with ingredients from N-body simulations a specified HOD 
makes strong predictions about a wide array of galaxy clus- 
tering statistics. The formalism thus allows us to use observa- 
tions of galaxy clustering to constrain the connection between 
galaxies and their host dark matter halos at each z, and in par- 
ticular to show that the luminous galaxies in the NDWFS have 
undergone significant merging or disruption between z — 0.9 
and z — 0.5 by comparing the HOD inferred from the z — 0.5 
clustering data to one which is passively evolved from the in- 
ferred HOD at z ~ 0.9. 

3.1. Simulations and mock catalogs 

Our modeling of galaxy clustering is based on mock cat- 
alogs constructed within the HOD framework by populating 
halos in a cosmological N-body simulation. We use a high 
resolution simulation of a ACDM cosmology (£Im = 0.25 = 
1-0 A , V B = 0.043, h = 0.72, n s = 0.97 and cr 8 = 0.8). The 
linear theory power spectrum was computed by evolution 
of the coupled Einstein, fluid and Boltzmann e quations us- 
ing the code described in IWhite & Scottl d!995l). This code 
agre es well with CMBf ast dSeljak & Zaldarriagall 19961) . see 
e.g. ISeljak et aH (120031) . The simulation employed 1024 3 
particles of mass 8 x NPIi^Mq i n a periodic cube of side 
500/i _1 Mpc using a TreePM code dWhitdl2002l) . The Plum- 
mer equivalent softening was 18/i _1 kpc (comoving). 

For each output we generate a c atalog of halos us ing the 
Friends-of -Friends (FoF) algorithm dDavis et al.f T985) with a 
linking length of 0. 168 x the mean inter-particle spacing. This 
procedure partitions the particles into equivalence classes, by 
linking together all particle pairs separated by less than a 
distance b. The halos correspond roughly to particles with 
p > 3/(27rZ? 3 ) ~ 100 times the background density. Our mass 
definition uses the sum of the particle masses in the halo, 
however to obtain better correspondence between our defi- 



nition of ha lo mass and that implic itly defined by the mass 
functions of lSheth & Tormenl d 19991) and lJenkins etaT] d2000t) 
we rescaled the masses by M/M fof = 1+0.01 (lnM fo f-23.5) 
where Mm is the FoF mass in units of /t'Mq . With this redef- 
in ition the mass function i n th e simulation lies bet ween those 
of lSheth & Tormenl d 1999!) and lJenkins etail d2000l) . differing 
from them by less than 10% in the mass range of interest. 

To make mock catalogs we use a halo model which distin- 
guishes between central and satellite galaxies. We choose a 
mean occupancy of halos: N(M) = (A'gai(Mhaio))- Each halo 
either hosts a central galaxy or does not, while the number 
of satellites is Poisson distributed about a mean A^at- With 
the luminosity-threshold sample s, we parameterize N(M) = 
N cen +N sat with 4 parameters (e.g. Zheng et al. 2 0051) 



N C en(M) = - erfc 



ln(M cut /M) 
V2a 



and 



AUM) = 



M-M a 
Mi 



(1) 



(2) 



for M > M cut and zero otherwise. Different functional forms 
have been proposed in the literature, but the current form 
is flexible enough for our purposes. Including a different 
low m ass r oll-off in the satellite term, fol l owing iTinker et alJ 
d2005l) and lConroy. Wechsler & Kravtsovl d2006h . does not al- 
ter our basic conclusions. 

Given an HOD and the halo catalogs we can produce a 
mock catalog in one of two ways. Central galaxies always live 
at the minimum of the halo potential . We either place satellite 
galax ies assuming an NFW profile dNavarro. Frenk & White! 
1996) with a concentration-mass relation fit to the halos in 
the simulation or anoint n sat dark matter particles, chosen at 
random, as galaxies. The two methods produce very simi- 
lar, though not identical, clustering, with the biggest differ- 
ences on Mpc scales. An analytic model (described in IZhengl 
l2004t ITinker et all 12005) also produces very similar results. 
The differences between the methods are smaller than the ob- 
servational errors, so we shall neglect them henceforth. 

3.2. Comparing with data 

From the model galaxy positions we compute £(r) in real 
space by direct pair counting in the periodic box for separa- 
tions < 20/T 1 Mpc. Beyond 20/r'Mpc we extrapolate assum- 
ing a constant bias. The redshif t distribution i s used to convert 
£(r) into w(9) using Eq. (50) of|Simon|(|2006), giving the pre- 
dicted clustering for any set of HOD parameters. We fit to the 
data assuming Gauss ian errors with the covariance matrices of 
iBrown et alJ d2006b). and assume a 5% error on the number 
density of galaxies. Figure [2] compares the best fitting HOD 
model predictions to the data at z — 0.9 and 0.5. 

In order to propagate the observational errors into uncer- 
tainties in the HOD (Eas. lll2b we used a Markov chain Monte- 
Carlo method (e.g. [G ilks, Richard son & Spiegelhalterlll996l) 
as detailed in IBrown et alJ d2006bl) . We found that the data 
were unable to rule out models with a 1 and a <C 1 , which 
we regard as unlikely for large red galaxies, so we impose a 
prior which penalizes a > 1 and a ~ 0. Tests indicate that sur- 
veys twice as large would not need this prior, though with this 
prior even our chains converged well. Because the mock cat- 
alog generation using NFW profiles is very fast and requires 
little memory we use this to generate the chains. 
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FIG. 2. — The angular correlation function, w(9), for the 0.4 < z < 0.6 and 
0.8 < z < 10 slices. Open circles with error bars represent the Bootes data, 
the solid line is the best fitting HOD model prediction and (in the lower panel) 
the dashed line is the prediction of the best fitting z ~ 0.9 model passively 
evolved to z — 0.5. 

3.3. Passive evolution 

For a subset of 250 of the models at z — 0.9 we use the 
particle-based method to produce mock catalogs which we 
passively evolve to z — 0.5 simply by tracking the particles 
based on their ID. The positions and halo memberships of 
these tagged particles are followed. We expect the HOD to 
be more robust than the small-scale clustering for these par- 
ticles 2 , but we show the latter in Figure |2] for completeness. 
Looking at the difference in the HODs is also more informa- 
tive, and gives us a clue as to what physics may be missing 
from pure passive evolution. 

4. DISCUSSION 

A comparison of the HOD of the passively evolved sam- 
ples with the HODs which best fit the z — 0.5 data indicates 
that evolution produces too many galaxies in high mass halos, 
as shown in Figure [3] A similar conclusion can be reached 
by comparing the clustering of the passively evolved mod- 
els to the data in Figure [2] — the excess clustering on small 
scales from passive models clearly indicates that the models 

2 The small-scale clustering depends on the evolution of the subhalos in- 
side of the host halo, and due to finite force and mass resolution these may 
not be correctly modeled in massive halos. 
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FIG. 3. — HODs for the z — 0.5 sample. The shaded area indicates the 
mean and standard deviation in the HOD from the Markov chains, fit to the 
w(6) data at z ~ 0.5. The mass scale has been increased by 7% as described 
in the text. Open triangles indicate the HOD from the models that fit the 
z ~ 0.9 data evolved to z — 0.5 by tracking particles. Squares assume 33% of 
the satellites in the passively evolved mocks disappear between z, ~ 0.9 and 
z~0.5. 

overpredict galaxy pairs within the same halo, i.e., they pre- 
dict too many satellite galaxies in high mass halos. Another 
indication of the excess is that the satellite fraction in the pas- 
sively evolved models is 0.24 ± 0.02 while that in the best 
fitting models is 0.18 ±0.02. 

There must be some physical process which reduces the 
number of galaxies in massive halos, and the most natural 
candidates are dynamical friction and tidal stripping which 
act to merge massive satellites with the central galaxy or dis- 
rupt them. Based on Figure [3] at high halo masses ~ 1/3 of 
the satellites in the passively evolved catalogs must have dis- 
appeared by z — 0.5. We caution that our calculation ignores 
sources and sinks. If galaxies can enter or leave the sample 
due to rapid evolution of their star formation rates, interpret- 
ting the evolution of galaxy clustering is non-trivial. However, 
z < 1 blue galaxies with masses comparable to the most mas- 
sive red galaxies are rare, and (at least at z — 0) red galaxies 
have little cold gas to fuel renewed bursts of star formation. 

There is a subtlety to bear in mind with these statistics. 
Merging results in a small fraction of z — 0.9 galaxies disap- 
pearing by z — 0.5. So we should really compare the passively 
evolved catalog (with n = 10~ 3 /i 3 Mpc~ 3 ) to a fitted one with a 
slightly lower number density - the "true" z — 0.5 descendants 
of the z — 0.9 galaxies. The cut-off mass scale for this catalog 
would be slightly larger than our fits and we should shift the 
mass definition in the HOD accordingly. We estimate the size 
of this shift, by considering the evolution of the central galax- 
ies, to be « 7% and in Figure[3]we increase the mass scale of 
the HOD inferred from the z — 0.5 data accordingly (though 
this shift does not change our conclusions). 

To estimate the impact of the disappeared satellites on the 
growth of the central galaxy or the boost of the intracluster 
light (ICL) in the halo, we make the following simple model. 
By matching the number density of halos above mass M with 
that of galaxies abo ve luminosity L (we use the B-band LFs of 
iBrown et alJ2006ah . we can relate the central galaxy luminos- 
ity to halo mass (the lumino us end of the LF is dominated by 
central galaxies as shown in Zheng et al. 2005). For our sam- 
ple, at z — 0.5 a halo of 5 x 1O 13 /i _1 M on average hosts a 
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central galaxy of sa 5L*, with L cen cx M for massive halos. 
If we also assume that the satellite LF has the same shape as 
the global LF, we can use the fitted HODs to find the number 
of satellites and integrate to measure their total light. We find 
that in such halos satellites on average have a total lumnosity 
of 2.6 L*, so if this is what is left after 1/3 of the satellites dis- 
appeared, the satellites would have contributed « 25% of the 
current stellar mass to the central galaxy or a similar mass to 
the ICL. For halos of 10 14 /i -I M Q the contribution is w 40%. 
These numbers are upper bounds, since satellites may have 
lost just enough mass to leave our sample. Note that the stel- 
lar material resides in the inner regions of the halo, where 
the potential well is the deepest, and is the last material to 
be disrupted. An "average" satellite would need to lose 40% 
of its stellar mass, or decrease its surface brightness, to drop 
out of the sample. Analyzing samples with lower luminosities 
would help us constrain this. 

We are unable, without further modeling, to differentiate 
between satellites fueling the growth of the central galaxy 
or the ICL, which at lower redshift co mprises ~ 5 - 10% of 
the stellar mass in groups and clusters dGonzalez et al.ll2.Q0Ql: 
IZibetti et all 120051: lAguerri et al.ll2006l) . We can, however, 
form an L cen -M relation at z — 0.9 as above. By matching 
progenitor halos to their descendants at z — 0.5 and assuming 
0.48 fi-band magnitudes of fading we find w 10% growth in 
the stellar mass of the central galaxy between z — 0.9 and 0.5. 
This would suggest the satellites also build an ICL compo- 
nent (the total stellar mass in the disappeared satellites, if it 
all ended up as ICL, constitute ~ 15% of the total halo lumi- 
nosity above the threshold). However this result relies on our 
assumption of a uniform LF and on stellar population evolu- 
tion predictions, and needs to be constrained by observations 
of the ICL. If we argue that there is little extended light at 
high z, either the satellites may have lost just enough mass to 
leave the sample or the "extra" mass would accrete onto the 
central galaxy. We regard this dichotomy as an open question 
requiring further investigation. 

Our conclusions are necessarily tentative due to the limited 
volume of the NDWFS Bootes survey, which does not probe 
the mass function above 10 14 /i _I M Q well. There are several 

3 The wider area survey need not be as deep as NDWFS. 



areas where more or different data would be beneficial. Tests 
using models of larger surveys 3 indicate that doubling the sur- 
vey volume removes the islands of parameter space which we 
have excluded with priors and shrinks the errors on the HOD 
parameters by « \f2. A measurement of the space density of 
groups richer than several members would shrink the errors 
on the high mass end of the HOD dramatically, but would 
require more volume than we have at present to contain a rep- 
resentative sample of rich groups. We also investigated the 
dependence of our results upon cosmology using similar sim- 
ulations with different parameters. Our results remain robust 
within the currently allowed range of models. 

This preliminary investigation shows the power of cluster- 
ing measures to inform questions of the formation and evo- 
lution of galaxies. We find evidence for evolution in the red 
galaxy HOD very different than that predicted by pure passive 
evolution models. Our result is largely independent of models 
of red galaxy stellar populations and estimates of dynamical 
friction scales. With more data from the NDWFS and future 
surveys we hope to be able to trace in detail the formation 
history of the most massive galaxies in the Universe. 
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